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Abstract 

Exact solutions for a model with variable G, A and bulk viscosity are obtained. Inflationary solu- 
tions with constant (de Sitter-type) and variable energy density are found. An expanding anisotropic 
universe is found to isotropize during its course of expansion but a static universe is not. The grav- 
itational constant is found to increase with time and the cosmological constant decreases with time 
as A oc t~ 2 . 
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1 Introduction 

In a recent paper, Kalligas, Wesson and Everitt (Kalligas et al 1995) have investigated a flat model 
with variable gravitational (G) and cosmological (A) "constants" . In the same line, Singh et al. (Singh 
el al 1998) have considered a viscous cosmological model with variable G and A. They considered 
a different energy conservation law from that of Arbab (Arbab 1997) but, they have found similar 
solutions as in (Arbab 1997). Very recently, we have studied Bianchi type I model and we have shown 
that the universe isotropize during its course of expansion (Arbab 1998). With a similar approach, 
we wish to study the effect of anisotropy in the universe with the energy conservation advocated by 
Singh et al. where G and A vary with time. We have shown that the introduction of bulk viscosity 
enriches this cosmology. 

Kalligas et al. have found solutions for a static universe with zero total energy density while G and 
A are allowed to vary with time. However, their solution, does not seem to be physically sensible, 
since one does expect G to vary with time in an empty universe ! 

In the present case, we show that the static universe must be empty and has a vanishing cosmological 
constant (A) and bulk viscosity (rf). We have also shown that the presence of viscosity helps an 
anisotropic universe to isotropize with expansion. We have also obtained solutions with constant and 
variable energy density that correspond to either static or inflationary universe. We remark that 
these solutions do not hold for our earlier work (Arbab 1998). The gravitational constant is found to 
increase with time. For a universe in balance (flat) expansion must accelerate in order to overcome 
future collapse. The presently observed acceleration of the universe is justified in the present models. 

2 SOLUTIONS FOR THE ISOTROPIC UNIVERSE 



In a flat Robertson Walker metric 

ds 2 = dt 2 - R 2 (t)(dr 2 + r 2 d6 2 + r 2 sm6 2 d(j> 2 ) (1) 
Einstein's field equations with time-dependent G and A read (Weinberg 1971) 

?V - ^g^n = 87iG(t)T^ + A(t)g^ (2) 

The variation of the gravitational constant was first suggested by Dirac (Dirac 1937) in an attempt 
to understand the appearance of very large numbers, when atomic and cosmic worlds are compared. 
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He postulated that the gravitational constant (G) decreases inversely with cosmic time. On the other 
hand, Einstein introduced the cosmological constant (A) to account for a stable static universe, as 
appeared to him at that time. When he later knew of the universe expansion he regretted its inclusion 
in his field equations. Now cosmologists believe that A is not identically but very close to zero. They 
relate this constant to the vacuum energy that first inflated our universe making it expanding. A 
vacuum energy could correspond to, from a point view of particle physics, a residue of some form 
of quantum field that is diluted to its presently small value. However, other cosmologists dictate a 
time variation of this constant in order to account for its present smallness. The variation of this 
constant could resolve some of the standard model problems. Like G the constant A is a gravity 
coupling and both should therefore be treated on an equal footing. A proper way in which G varies is 
incorporated in the Brans-Dicke theory (Brans & Dicke 1961) . In this theory G is related to a scalar 
field that shares the long range interaction with gravity. In the literature A takes several ansatz, like 
A cx R~ 2 , A tx H 2 , A tx 8ivGp, etc, and with different reasoning (Overdin & Cooperstock 1998 and 
Sahni & Starobinsky 1999). In the present work, we allow the variation of G to be cancelled by the 
variation of A (Beesham 1986, Abdel Rahman 1990, Pande 2000 Bonanno & Reuter 2002). 
Considering the imperfect-fluid energy momentum tensor 

Tpv = (p + p*)u^u v - p*g fll , (3) 

eq.(2) yields the two independent equations 

= -^G(3p* + p) + A , (4) 



R 



and 



7T 2 



3 ( - ) = SirGp + A . (5) 



Elimination of R and the differentiated form of eq.(5) gives 

A \ 

(6) 



3(p- + p)A=-(§p + p +i A,)* 



where a dot denotes differentiation with respect to time t and p* — p — Sr/H , rj being the coefficient of 
bulk viscosity, H the Hubble constant. The equation of state relates the pressure (p) and the energy 
density (p) of the cosmic fluid by the equation 

p=( 7 -l)p (7) 

where 7 = constant. Vanishing of the covariant divergence of the Einstein tensor in eq.(2) and the 
usual energy-momentum conservation relation (Tl£" = 0) lead to 

8ttG P + A = , (8) 

and 

P + 3( P * +p)H = 0, (9) 

or 

p + 3{p + p)H = 9 V H 2 . (10) 

We see that the bulk viscosity appears as a source term in the energy conservation equation. Hence 
the RHS of eq.(10) would correspond to the rate of thermal energy generated due to viscosity. This 
may help solve the generation of entropy in the universe associated with the standard model of 
cosmology. 

In this paper we will consider the very special ansatzs (Arbab 1997) 

A = 3f3H 2 , (3 = const. , (11) 
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and 

V — Vop" , Wo > 0, n = const. (12) 

We have shown very recently that eq.(ll) is equivalent to writing A as A = (^3) 4nGp for a non- 
viscous model (Arbab 2002). This form is interesting since it relates the vacuum energy directly to 
matter content in the universe. Hence, any change in p will immediately imply a change in A, i.e. if 
p varies with cosmic time then A also varies with cosmic time. 

In what follows we will discuss the solution of the model equations for an isotropic and an 
anisotropic universe with the above prescription for A and r). 

2.1 Solution with constant energy density 

One can satisfy eq.(10) with a constant energy density (p = const) with: 

(i) H — 0, which implies a static universe. 

(ii) 77 = rjop (i.e. n = 1), H = = const.. The solution of this equation is of the form R — 
const. exp(Ht). We remark here that the classical inflation with an equation of state p — —p is not 
permitted in this model. 

2.2 Solution with variable energy density 

Inflationary solution with variable energy density 

Consider the bulk viscosity to have the form 77 = r/op (i.e. n = 1). With H = Ho < ^ we have 
R — const. exp(Hot) so that eq.(10) yields a decaying mode of the energy density given by 

p = Fexp -3^0(7 - 3r) H )t , F = const.. (13) 

Now consider the case (3 = 1. 
Equations (5), (8) and (11) yield 

A = const. , G = (14) 

We remark that this solution is not possible within the framework of the conventional inflationary 
models. It is however remarked by Abdel Rahman (Abdel Rahman 1990) that one possible inter- 
pretation of eq.(14) is that the universe came to being just prior to the onset of gravity at t = as 
a result of a vacuum fluctuation propelled by the repulsive effect of the positive cosmological constant. 



3 SOLUTIONS FOR AN ANISOTROPIC UNIVERSE 

For the Bianchi type I metric 

ds 2 — dt 2 — R 2 dx 2 — R\dy 2 — R^dz 2 (15) 
with an imperfect-fluid energy momentum tensor, Einstein's field equations yield (Arbab 1998) 

R1R2 R1R3 R2R3 „ „ . , 1c s 
R^ + R^R 3 + Rjir- 8nGp - A ' (16) 

(17) 
(18) 
(19) 



R1R2 


+f 


+ R2 

R2 


= 8ttG P * 


-A 


R1R3 

R1R3 




+ fh 

R3 


= S-kGp* 


-A 


R2R3 

R2 R3 


4 


+ R3 

Rs 


= 8irGp* 


-A 
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and 



SivGp + A + 8nG 



, Ri R2 R3 



= . (20) 



where a dot denotes differentiation with respect to time t. From eqs.(16)-(20) one obtains 

^- + ^ + ^=^G{p + 3p')-A. (21) 

Hi H2 -K3 

where p* = p — 3r/H. The energy conservation (T^J = 0) implies that 

p + 3H(p + p*) = , (22) 

or 

p + 3H(p + p) -9riH 2 = . (23) 
Here we define the average scale factor R by R = {R\R2R 3 ) 1 ^' i so that 

*=! = U#+#+#Y m) 



R 3 \Ri R 2 Rs^ 
Using eqs.(22) and (24), eq.(20) yields 

A + 87rGp = . (25) 
Let us now assume that the energy density is given by the power law 

p = At m , A — const. , m= const. , (26) 

and the average scale factor 

R = Bt a , a = const. , B = const.. (27) 

Substituting eqs.(26) and (27) in (23), (11) and (25) one gets 

p= At- 1/(1 - n) , (28) 

A = 3af3t~' 2 , (29) 
3/3a 2 (l-n) (2n _ 1)/(1 _ n) 1 
°" 4jrA(2n-l)* ' ^2 ' 1 ' m 

and the condition m = j^. For physical significance G > so that n > |. This implies that 
the gravitational constant is an ever-increasing function of time. Consequently, for a flat (balanced) 
universe the expansion must increase (accelerate) so that the universe can remain in balance. Thus, 
the presently observed acceleration of the universe may be attributed to this ever growing gravity 
instead of invoking any exotic matter. 
We now consider the anisotropy energy (a) defined by 

Using eqs.(5), (16) and (24), the above equation becomes 

87rGa = 18^"j +48ttGp + 6A. (32) 



We see that the anisotropy energy (a) becomes 

8vrGa = Dr 

a result that has been obtained by Arbab (Arbab 1998) 



87rGcr = Dt 2 , D = const. , (33) 
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3.1 Constant energy solution 

Using eq.(12), eq.(23) can be written in the form 

p + 3-yHp = 9ri oP n H 2 . (34) 

Now consider the following two cases: 

(i) Static universe ( H = , R = const.). 

The above equation yields p — const, and eqs. (11) and (25) yield A = and G = const.. It follows 
from eqs. (5) and (32) that the anisotropy energy a = 0. It is shown by Kalligas et al. ( Kalligas et 
al. 1995) that a static universe can only be isotropic, i.e., a — . However, a static universe with 
constant energy density can not exist unless G — 0, as is evident from eq.(5). Hence, only an empty 
static universe can exist. 

We remark here that the claim made by Kalligas et al. that a universe with a vanishing total energy 
density has G and A vary with time is physically nonsensical. 

(ii) Inflationary universe with constant energy density (p — const. , H = const., n = 1) 
H = Ho — const, implies R — const. exp(Hot). Equation (34) gives Ho = ^- . 

Hence, eqs. (11), (25) and (32) give A = const., G = const, and a — (j^^J = ^tt^g ) ' wnetner 
A = or not. 

4 Conclusions 

In this paper we have studied both isotropic and anisotropic models with variable G, A and bulk 
viscosity. We have found that energy conservation is guaranteed provided the three scalars, G, A and 
7] conspire to satisfy it. We have found that a constant energy density would lead to either a static 
or inflationary universe. We have found that the classical inflation of an equation of state p — —p 
is not permitted. Inflationary solutions with constant and variable energy density are found. These 
solutions are influenced by the presence of the bulk viscosity. Such a solution is triggered by the 
presence of the bulk viscosity alone. An initially expanding anisotropic universe is found to isotropize 
as it evolves. However, a static empty universe must be isotropic. This may explain why the present 
universe is very isotropic. Finally, we have found that the gravitational constant to increase with 
time in radiation-dominated (7 = |) and matter-dominated (7 = 1) epochs. The presently observed 
acceleration of the universe may be interpreted as due to a balance between expansion and gravity. 
And since gravity increases expansion must accelerate! 
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